Critical scaling in standard biased random walks 
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The spatial coverage produced by a single discrete-time random walk, with asymmetric jump 
probability p 1/2 and non- uniform steps, moving on an infinite one-dimensional lattice is investi- 
gated. Analytical calculations are complemented with Monte Carlo simulations. We show that, for 
appropriate step sizes, the model displays a critical phenomenon, at p = pc- Its scaling properties 
as well as the main features of the fragmented coverage occurring in the vicinity of the critical point 
are shown. In particular, in the limit p pc, the distribution of fragment lengths is scale-free, with 
nontrivial exponents. Moreover, the spatial distribution of cracks (unvisited sites) defines a fractal 
set over the spanned interval. Thus, from the perspective of the covered territory, a very rich critical 
phenomenology is revealed in a simple one-dimensional standard model. 
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Since the beginning of the past century, random walk 
(RW) theory has allowed to deal with a diversity of prob- 
lems in a number of areas of physics, as well as in many 
other theoretical and applied fields, e.g, biology, chem- 
istry, computer sciences and finance [l| . The undoubtful 
importance of RW models, with their wide range of dis- 
tinct applications, stems from their simplicity and effec- 
tiveness in modeling systems experiencing disorder, noise 
or randomness, which are ubiquitous features of real sys- 
tems. In particular, in physics, RWs can be seen as the 
"harmonic oscillator" of disordered and stochastic sys- 
tems, serving as starting point for more realistic models. 

A fundamental quantity in any phenomenon where 
RWs are relevant is the number of distinct sites visited, 
since it furnishes the extent of the active territory. In- 
deed, it is crucial in processes ranging from reaction ki- 
netics to population dynamics, and also in technical ap- 
plications such as in search strategies [l, Q . As a conse- 
quence, analytical and numerical estimates of the covered 
territory are available for lattices of different geometry, 
dimensionality and boundary conditions 0, 0] ^ for diverse 
statistics of jumps, symmetric or not J6|, and other vari- 
ants • Time covering pro blems [1, [J and coverage by 
a large number of RWs [l0| have been investigated too. 
The vast literature on coverage mainly deals with two 
dimensions, although there are also many works about 
the standard symmetric one-dimensional (ID) RW (e.g., 
[IjHQ)- Meanwhile, as far as we know, little or no atten- 
tion has been paid to the asymmetric ID case, despite of 
its importance in biased or anisotropic processes such as 
electrophoresis, polymer translocation through pores and 
Brownian ratchets. However, as we will show, the asym- 
metric ID problem presents its own peculiar features and 
nontrivial scaling properties. 

In the present work, we investigate the coverage of an 
infinite ID regular lattice by a single RW characterized 
by: i) asymmetry, that is, at each independent step there 



is a probability p ^ 1/2 to step, let us say, to the right, 
and additionally, ii) distinct step sizes in opposite direc- 
tions. Let us call Z+ and l~ the sizes of the steps in the 
positive and negative directions, respectively. They will 
be expressed as integer multiples of the arbitrary lattice 
parameter. In the symmetric case = l~ , only the po- 
sitions that are multiple of l'^ are reachable. Moreover, 
the covered fraction of the interval spanned by the RW 
is l//"*", independently of p. In particular, if Z+ = Z~ = 1, 
complete coverage of the RW span occurs. However, for 
the asymmetric case ^ l~ , where the two anisotropic 
ingredients compete, a nontrivial changeover between dif- 
ferent coverage regimes, dependent on p, may take place. 
In fact, we will show that a critical phenomenon occurs 
as the jump probability p reaches a critical value. More- 
over, we will characterize the transition as well as the 
partially covered, fragmented, states, focusing on their 
scaling properties. 

The general basic outlines to determine the number 
of distinct sites visited by a RW can be found, for 
instance, in Refs. In general, the average number 

of different sites visited at step n, S'„, can be expressed 
'AS Sn = 1 + Ss^oSr^i Pii^)' where F^(s) is the 
probability that the walker arrives at site s for the 
first time at step i. Moreover, Fi{s) and Pj{s) (the 
probability that, at time step j, the walker is located at 
integer position s) are related through 
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F,{s)Pn-^{Q), for7i>l. 



(1) 



Pois) = Ssfi- Then, from Eq. ([T]), one ob- 
the following relation between generating func- 
P(s,z) = S,fi + F{s,z)P{0,z), where 



while 
tains 
tions: 

Pi^,^) = En>oPn{s)z^ and F{s,z) = E„>l^^nW^ 
Assuming \z\ < 1, one obtains 



Siz) ^ [{l-zfPiO,z)y 



(2) 



where S{z) = J2n>o 



For the present problem, it is 
easy to show that P(0, z) explicitly is 
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I- i+ 
with z = zp^^+'^ (1 — p) . From the definition of 

S{z), the quantity Sn can be obtained as 1/n! times the 
nth derivative of S{z), evaluated at z = 0. 

If and l~ have common factors, a mapping ex- 
ists into the corresponding case of reduced (mutually 
prime) lengths. Therefore, we will restrict our study to 
asymmetric coprime couples of step lengths. Within the 
latter class of RWs, one has the subclass where one of 
the lengths is unitary. Let us consider as representative 
of this subclass, the case {l^,l~)= (2,1) that admits an 
exact solution. In this case, the sum in Eq. jS]) becomes 

^(0, z) - Ek>o ( ? ) /(I - P)^'^'^ that can be 
reduced to 

P(0,z) = K(z2;+v/T^)i/VyT^, (4) 
for \y\ < 1, where = (27/4)p(l — p^z^. Taubcrian 
methods can be applied to evaluate Sn 0]. Alterna- 
tively, the nth derivative of S{z) can be calculated 
through Cauchy integral formula over a suitable 
contour encircling the origin. Since S{z) given by 
Eq. ^ has one single pole in the complex plane, 
at z = 1, then, in the limit of large n, one gets 
(after conveniently deforming the integration path) 

Sn = ~d [2"+ip(0, z)] /dzU=i = (n + 1)/P(0, 1) + Co, 
where cq is a constant of order 1. It is noteworthy that 
this is the same asymptotic law found for the stan- 
dard RW, with unbiased symmetric jumps to nearest 
neighbors (hence l'^=l~=l), but in 3D regular lat- 
tices Q . The fraction of different sites visited (measured 
over the average length of the RW) is /u^„ = S'„/L„, 
where Ln is the average total displacement. In the 
large n limit, the length of the RW, for p ^ 1/3, is 
Ln ~ |(s)n| = |3p - l|n. Thus, asymptotically, /i,,„ 
becomes /„ = [\3p — 1| P{0, hence, the fraction of 
unvisitcd sites is 

fu = l-fv = l-[|3p-l|P(0,l)]-\ (5) 

where P(0, 1) is given by Eq. 

Fig. [Ija) exhibits /„ as a function of p, for (Z+,Z^)= 
(2,1). A transition occurs at pc = 1/3, where /„ van- 
ishes as /„ = A -I- ©(A^), with A = p — pc, that can 
be derived exactly from Eq. ([5]). For p < Pc all sites are 
eventually visited at least once, as expected, because, as 
soon as (s)„ = 3A < 0, the walker is biased towards 
the direction of unitary steps, which in turn implies full 
coverage of the RW length. Meanwhile, for p > pc, se- 
quences of adjacent visited sites (fragments) arc inter- 
rupted by unvisitcd ones. Therefore, the RW undergoes 
a transition from a fully covered state to a fragmented 
one. For other instances of (Z+,1), the transition oc- 
curs at the critical probability pc ~ l/(^^ + 1); where 
{s)n = {pl~^ + P — changes sign (driftless diffusion). 
The case (Z+,^^)= (3,1), obtained by means of Monte 
Carlo (MC) simulations up to n ?» lO'' time steps, is also 
displayed in Fig.[IJa), exhibiting similar features. In both 
cases, /«(A) vanishes with unitary exponent (see inset of 
Fig. Ha)). 
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FIG. 1: Fraction fu of sites left unvisited as a function of p. In 
all cases, symbols correspond to MC simulations, and dotted 
lines are guides to the eye. (a): {l'^,l~)= (2,1) (circles) and 
(3,1) (squares). The full line corresponds to the theoretical 
prediction given by Eq. ((SJ. Inset: fu vs. A = p — pc in log- 
log scale for the same data of the main frame, (b): (/"*",/") 
takes diverse coprime values indicated on the figure. 

For non- unitary coprime step lengths (see Fig. [Hb)) 
a more general scenario arises. Full coverage occurs 
only at the critical point Pc = + l~), where 

{s)n = {pl^ + {p — ^)l~)n is strictly null. Fragmented 
states are found both below and above Pc, with maximal 
unvisited fractions, = 1 — l/?~ and /+ = 1 — 
respectively. Thus, the cases (^^^,1), with Z+ > 1, con- 
stitute special instances where one of the states is fully 
covered, in accordance with the fact that the correspond- 
ing maximal unvisited fraction /" vanishes. Although we 
are not dealing with symmetric steps, notice that in the 
symmetric case (1,1), /~ = = and the full curve 
fu{p) collapses to zero, in agreement with the facts that 
there is no transition in such case and that full coverage 
occurs for any p. 

As a paradigmatic example, we will analyze the 
analytically soluble case (2,1), in the vicinity 

of the critical point, i.e., in the limit A ^ 0+. In order 
to quantitatively characterize fragment sizes, the usual 
computed quantities are [ll] |: 

m = ^^i' = E "^i^y (6) 
e>i e>i e>i 

where ni is the mean number of fragments of size £, 
normalized per site. Since two contiguous fragments are 
separated, in the (2,1) case, by one single unvisited site, 
then hi « /„, that vanishes as ~ A (sec Fig.[TJa)). Also, 
straightforwardly, X]£>i ~ ^ ~ fui that approaches 1 
in the critical limit. Noticing that ng£ is the probability 
that a given site belongs to a fragment of size £, then, 
£ w S^>i'^^^^ defines the mean size of the fragments. 
In order to compute {£), the distribution of sizes of 
covered clusters (or fragments), n^, was numerically 
built from MC simulations run up to n sa 10^/ A steps 
and averaged over at least 10^ different realizations. 
The distributions for different values of A are displayed 
in Fig. [21 For very large £, the decay is exponential: 
^ exp(--^/A). Parameter A, together with {£), are 
plotted as a function of A in the upper inset of Fig. [2] 
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(being A w (£)/2 ~ A""^, with 7 « 1.15). Meanwhile, 
rii ~ A, representing a finite fraction of /„. In the lower 
inset of Fig. [21 the same distributions of the main frame 
are scaled. Let us employ the standard ansatz for cluster 
size distributions [ll|, defined through: 

nf(A) oc A'^0(Ai/'^£)/(Ai/'^^)^, (7) 

where (j){x) goes to a constant value for small x and de- 
cays exponentially in the opposite limit of large x. The 
power-law decay, with exponent r « 1.15, that emerges 
in the limit of vanishing A is characteristic of a critical 
behavior and signals the coexistence of fragments of all 
sizes in that limit. 




10' 10^ 10^ 10* „ 10= 



A^~"j', with ttp = 1, as in the present problem. On the 
other hand, (£) A^'^'p, with 7p = 1 and Wp = = 2, 
values that are close but different from those found for the 
present problem. Moreover, the distribution of fragment 
sizes is a power-law, in contrast with the pure Poissonian 
one for IDP. Then, we may conclude that the present 
model docs not belong to the IDF universality class. 
Indeed, by identifying visited sites with occupied ones, 
the occupation probability in our problem is /i,, that 
tends to one in the critical limit. However, differently 
from the standard percolation problem, in the present 
case, unvisited sites are not independently located, e.g., if 
^~) = (2,1), a sequence of two or more adjacent unvis- 
ited sites has associated a strictly null probability of oc- 
currence. Therefore, occupation correlations arise which 
are absent in the standard percolation problem. 

Concerning unvisited sites, their spatial distribution 
was investigated through a box-counting procedure [TTI . 
[T^ . From the history of a single RW, a segment of length 
L = 2^" K, 10^ was divided into boxes of length 2^ , with 
fc > 0. For each e = 2'^/L, the number of boxes con- 
taining unvisited sites, N{e), was computed. Outcomes, 
accumulated over 10^ realizations, are displayed in Fig.O 
The neat behavior N{£) ~ e~'^f , for small e, means that 
the spatial distribution of unvisited sites constitutes a 
fractal set, with dimension df. Moreover, the fractal ex- 
ponent is in good accord with the exact scaling relation 

df=T-l. 



FIG. 2: Distribution of the sizes of covered fragments (ni was 
omitted), for {l^ ,1~) — {2,1) and different values of A = 10~^, 
3 X 10~^, . . ., 10~*, from a to g, respectively. Upper inset: 
mean size of fragments {I) (squares) , inverse exponential rate 
A (circles), and [lOni]^^ (triangles) as a function of A. Lower 
inset: Scaling plot of all the distributions represented in the 
main frame, with uj — 2/a and 1/a = 1.15 ± 0.05. Dashed 
lines are drawn for comparison and their slopes indicated on 
the figure. 

By means of integral approximations to the sums in 
Eqs. © and employing Eq. ([7]), one gets the follow- 
ing relations amongst critical exponents. Firstly, 1 « 
E£>i«f^ ~ ir^f^dl ^ A— 2/-, implying oj = 2/a. 
Secondly, A"^ ~ (i) « nei'^di ~ A'^-^/'", hence 

= 3/(7 — 7, that, together with the preceding relation, 
implies 7 = l/tr and cu = 2/a. The latter equality is in 
good accord with the behavior of the envelope of the dis- 
tributions that has slope -2 (Fig.[2|). Excellent data col- 
lapse is obtained for cu = 2/a, with 7 « 1.15. Addition- 
ally, since he ~ A, then, from he w nedl ^ A"^'^''', 
it must bcr = 2 — a = 2 — I/7. From the scaled his- 
tograms, we obtained r « 1.15 ± 0.05, consistent with 
the theoretical prediction within error bars. 

At this point, it is worth comparing our results with 
those for another ID critical phenomenon, namely ID 
percolation (IDP) with bonds connecting nearest neigh- 
bors [T^ . to which many important ID models are re- 
lated (e.g., Ref. [l^l)- On one hand, for IDP, hi ~ 
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FIG. 3: Scaling plot of the number of boxes A'^ containing 
unvisited sites as a function of e (box size in units of L, where 
L = 2^°) for (/+,r) = (2, 1) and different values of A indi- 
cated on the figure. The scaling exponent is(5 = l/(l — d/), 
where df ~ 0.15. Inset: original plots of the data scaled in 
the main frame. All solid lines are drawn for comparison and 
their slopes indicated on the figure. 

In summary, we have investigated the spatial cover- 
age of single discrete-time anisotropic RWs, moving on 
an infinite one-dimensional lattice. Anisotropy manifests 
both in the length (/+, l^), as well as in the probabilities 
(p, 1— p), of jumps in opposite directions. We revealed the 
existence of a critical phenomenon in ID that may result 
from the competition between the opposite trends pro- 
vided by the two anisotropic ingredients (step size and 
step probability). We illustrated our findings with the 
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FIG. 4: Critical behavior for (l^ ,l~) = {5,3). Results corre- 
spond to the limit A ^ 0"*", but the same exponents are found 
in the limit A — > 0~. (a) Scaling plot of the distribution of 
the sizes of covered fragments, values of A as in Fig. ((2)1, with 
Lo = 2/cr and 1/cr = 1.5 ± 0.2, r ~ 1.4. (b) Scaling plot of the 
number of boxes A*' containing unvisited sites as in Fig. ([Sjl. 
In this case d/ ~ 0.4. The dotted line in (a) and solid lines 
in (b) are drawn for comparison and their slopes indicated on 
the figure. 

particular case in which the steps are = (2,1), 

which undergoes a transition from fully to partially cov- 
ered states as the jump probability p overcomes a criti- 
cal value. The power-law distribution of sizes of covered 
segments, occurring in the limit p — > p^f, indicates the 
coexistence of fragments of all lengths, with no charac- 
teristic length scale. Moreover, the spatial distribution of 
scission points (unvisited sites) determines a fractal set, 
in contrast with other models where the deposition of 
cracks has common statistics (e.g., ID percolation [T^j, 



scission model [15|)- It is pertinent remarking that akin 
features have been observed in one-dimensional reaction- 
diffusion [3, [13] J ^-state Potts spin flipping and frag- 
mentation dynamics [isj , although criticality is attained 
as time evolves and critical exponents are different. A 
possible connection remains to be investigated. 

Other asymmetric instances with steps (/^,1), whose 
critical curves are illustrated in Fig.llfa), display a qual- 
itatively similar picture to the case (2,1). Meanwhile, if 
both steps take non-unitary coprime values (Fig. [^b)), 
the same critical phenomenology is observed in both lim- 
its A ^ 0^. As a further example, scaling plots are 
also displayed, in Fig. IH for the case (^''', Z~)=(5,3) in 
the limit A ^ 0+. In general, critical exponents related 
to the fractal dimension are not universal but depend on 
the step lengths, since distinct site occupation correla- 
tions take place. 

One one hand, the asymmetric RW, seen from the 
present perspective, may bear interest per se because 
of the nontrivial criticality contained in a simple model. 
On the other hand, it may constitute a useful statistical 
paradigm for the formation of domains or fragments by 
a non-equilibrium process driven by biased signal propa- 
gation. Additionally, the current coverage problem may 
be potentially useful in technical applications, e.g., in 
search strategies such as for cache hit/miss ratio opti- 
mization 

Acknoledgements: We acknowledge Brazilian agen- 
cies CNPq and Faperj for partial financial support. 



[1] G.H. Weiss, Aspects and applications of the random walk, 
(North-Holland, New York, 1994); B.D. Hughes, Ran- 
dom Walks and Random Environments, Vol. 1: Random 
Walks (Clarendon Press, Oxford, 1996); P.G. Doyle, and 
J.L. Snell, Random Walks and Electric Networks (The 
Math. Ass. of America, Inc., Washington-DC, 1999). 

[2] G.M. Viswanathan, S.V. Buldyrev, S. Havlin, M.G.E. da 
Luz, E.P. Raposo, and H.E. Stanley, Nature (London) 
401, 911 (1999); E.P. Raposo, S.V. Buldyrev, M.G.E. da 
Luz, M. C. Santos, H.E. Stanley, and G.M. Viswanathan, 
Phys. Rev. Lett. 91, 240601 (2003). 

[3] D. Thiebaut, IEEE Trans. Comput. 38, 1012 (1989); I. 
Gluhovsky, D. Vengerov, and B. O'Krafka, ACM Trans. 
Comput. Syst. 25, Art. 2 (February 2007). 

[4] E.W. MontroU, and G.H. Weiss, J. Math. Phys. 6, 167 
(1965). 

[5] J. Vineyard, J. Math. Phys 4, 1191 (1963); H. Larralde, 
and G.H. Weiss, J. Phys. A: Math. Gen. 28, 5217 (1995). 

[6] M. Ferraro, and L. Zaninetti, Phys. Rev. E 64, 056107 
(2001). 

[7] J.E. Gillis, and G.H. Weiss, J. Math. Phys. 11, 1307 
(1970); P. Molinas-Mata, M.A. Munoz, D.O. Martmez, 
and A.L. Barabasi, Phys. Rev. E 54, 968 (1996); F. van- 
Wijland, and H.J. Hilhorst, J. Stat. Phys. 89, 119 (1997); 
F. vanWijland, S. Caser, and H.J. Hilhorst, J. Phys. A 
30, 507 (1997); S.-Y. Huang, X.-W. Zou, W.-B. Zhang, 
and Z.-Z. Jin, Phys. Rev. Lett. 88, 056102 (2002); E. Al- 
maas, R.V. Kulkarni, and D. Stroud, Phys. Rev. E 68, 



056105 (2003). 

[8] K.R. Coutinho, M.D. Coutinho, M.A.F. Gomes, and 

A.M. Nemirovsky, Phys. Rev. Lett. 72, 3745 (1994). 
[9] M.S. Nascimento, M.D. Coutinho, and C.S.O. Yokoi, 

Phys. Rev. E 63, 066125 (2001). 
[10] H. Larralde, P. Trunfio, S. Havlin, H.E. Stanley, and G.H. 

Weiss, Phys. Rev. A 45, 7128 (1992); iMd., Nature 355, 

423 (2002). 

[11] J. Feder, Fractals (Plenum Press, New York, 1988). 
[12] P.J. Reynolds, H.E. Stanley, and W. Klein, J. Phys. A 10, 

L203 (1977); W. Klein, H.E. Stanley, S. Redner, and P.J. 

Reynolds, J. Phys. A 11, L17 (1978); D. Stauffer, and 

A. Aharony, Introduction to percolation theory (Taylor & 

Francis, London, 1992). 
[13] K. Christensen, H. Flyvbjerg, and Z. Olami, Phys. Rev. 

Lett. 71, 2737 (1993). 
[14] T. Tel, Z. Naturforsch 43a, 1154 (1988). 
[15] E. Ben-Naim, and P.L. Krapivsky, Physica D 107, 156 

(1997); and references therein. 
[16] G. Manoj, and P. Ray, J. Phys. A: Math. Gen. 33, 5489 

(2000). 

[17] B. Derrida, V. Hakim, and V. Pasquier, Phys. Rev. Lett. 
75, 751 (1995). 

[18] Z. Cheng, and S. Redner, Phys. Rev. Lett. 60, 2450 
(1988). 



